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ABSTRACT 

We show, using global 3D grid-based hydro dynamical simulations, that Ultra Fast Outflows (UFOs) 
from Active Galactic Nuclei (AGN) result in considerable feedback of energy and momentum into the 
interstellar medium (ISM) of the host galaxy. The AGN wind interacts strongly with the inhomoge- 
neous, two-phase ISM consisting of dense clouds embedded in a tenuous hot hydrostatic medium. The 
outflow floods through the inter-cloud channels, sweeps up the hot ISM, and ablates and disperses 
the dense clouds. The momentum of the UFO is primarily transferred to the dense clouds via the 
ram pressure in the channel flow, and the wind-blown bubble evolves in the energy-driven regime. 
Any dependence on UFO opening angle disappears after the first interaction with obstructing clouds. 
On kpc scales, therefore, feedback by UFOs operates similarly to feedback by relativistic AGN jets. 
Negative feedback is significantly stronger if clouds are distributed spherically, rather than in a disc. 
In the latter case the turbulent backflow of the wind drives mass inflow toward the central black 
hole. Considering the common occurrence of UFOs in AGN, they are likely to be important in the 
cosmological feedback cycles of galaxy formation. 

Subject headings: galaxies: evolution — galaxies: formation — hydrodynamics — ISM: jets and 
outflows — methods: numerical 



1. INTRODUCTION 

Active Galactic Nuclei (AGN) winds have for a long 
time been considered an integral part of the feedback cy - 
cle of galaxy formation flSilk fe Rees|1998||Fabian|1 999). 
The kinetic energy fed back by the wind into the inter- 
stellar medium of the host galaxy is, in principle, suf- 
ficient to heat, disperse, and possibly unbind dense gas 
and therefore inhibit galaxy-wide star formation ([Cren- 
shaw fc Kraemer|2012|) . Such a mechanism, whereby the 
central supermassive black hole (SMBH) is capable of 
controlling the growth of its host galaxy on scales much 
greater than its gravitational sphere of influence is an 
attractive idea to explain the apparent co-evolution of 
SMBH and galaxy as evidenced by the M-a re lation 
(iFerrarese & Merritt||2000| IGebhardt et al.||2000|), and 



the shape and evolution of galaxy lumino sity functions 
and BH mass functions over cosmic time (Croton et al. 



2006] |Croom et alT][2QQ9] |Merloni fc HeTnz] |2QQ8 ) . 



Mass outflows are common in all types of AGN and, 
those most closely associated with dis c winds are ob- 
serve d in UV and X-ray absorption lines (Crenshaw et al. 
2003). An extreme class or component of outflows re- 
cently detected in highly ionized and highly blueshifted 
Fe K-shell absorption lines in the hard X-ray band are 
ultra f ast outflows (UFOs, Cappi|[2"006| |Tombesi et al. 
2010b). UFOs are thought to be mildly relativistic disc 
winds, with speeds v ~ 0.01c to 0.1c (several 10 3 kms _1 
to 10 4 kms -1 ), originating within 100 - 10 4 gravitational 
radii o f the SMBH. From a sample of 42 local radio-quiet 
AGN, |Tombesi et"aLl Q2010a| |2012| ) determined that the 
incidence of UFOs in AGIN is greater than 40%, that mass 
outflow rates are typically 0.01 - 1 M y _1 , and that the 
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outflow kinetic power is 10 42 - 10 45 ergs _1 . These au- 
thors also suggested that such outflows ought to have a 
strong feedback effect on the evolution of the host galaxy. 

The effectiveness of any mode of AGN feedback de- 
pends sensitively on the properties of the ISM, but there 
are few detailed st udies incorporating r ealistic multi- 
phase distribu tions ( Saxton et al.[2005| [Sutherland fc 
| Bicknell| [2007) SB Thenceforth; JCooper et ai"H 2008 1 
Gai bler et al.||2012| |Wagner et aT1|2012[ WBU12 h ence- 
forth). WBU12 showed with 3D hydrodynamical simu- 
lations of relativistic AGN jets interacting with a two- 
phase ISM that jets with powers 10 43 - 10 46 ergs _1 are 
capable of efficient energy and momentum transfer to 
disperse the dense gas in the bulge of galaxies to veloci- 
ties commonly observed in (high-redshift) radio galaxies 
flMorganti et al |2005[ [Nesvadba et al.|2 008), and beyond 
those predicted by the M-a relation, if the Eddington 
ratio of the jets is greater than 10 -4 . The dominant 
force responsible for the efficient energy and momentum 
transfer was identified as the ram pressure carried by jet 
streams that percolate the porous two-phase ISM. 

Although UFOs, on average, have kinetic powers an or- 
der of magnitude less than AGN jets, their mass outflow 
rates are comparable to their accretion rates, and thus, 
for the same kinetic power, they c arry considerably more 
momentum ( |Tombesi et al.||2012| ) than jets. UFOs may 
therefore substantially affect the galaxy-scale ISM of the 
host, in particular the dense, warm and cold phases of 
the ISM from which stars could form. The present let- 
ter tests this proposition with global 3D hydrodynamic 
simulations of UFOs interacting on kpc scales with the 
two-phase ISM of the host galaxy. 

2. EQUATIONS, CODE, AND INITIAL AND BOUNDARY 
CONDITIONS 

Let p, v, p, I, T, A, 0, 7 = 5/3, /c, and \i be the 
fluid density, the three dimensional velocity vector, the 
pressure, the unit tensor, the temperature, the cooling 
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rate, the gravitational potential, the adiabatic index for 
an ideal gas, Boltzmann's constant, and the mean mass 
per particle, respectively. The system of equations de- 
scribing the UFO outflow, hot atmosphe re, and warm 
clouds in t he one-fluid approximation is ( |Landau fc Lif- 
shitz|l987| |: 
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= VA(T)-pV0.v; 
p = pkT/p. 

We integrate these equations using the publicly avail- 
able, open-source Eu lerian Godunov-type code PLUTO 
flMignone et al.||2007|) version 3.1.1. 

The UFO inlet is placed at x = (0, 0, 0) and is directed 
in the positive x-direction with an opening angle = 30° . 
The velocity, mass injection rate, and energy injection 
rate, at the inlet are v = 0.1c, Mufo = 0.1 M y -1 , 
and Pufo = 10 44 ergs -1 , respectively. These are typical 
UFO parameter values, and represent a case where the 
internal and kinetic energy injection rates are compara- 
ble. The energy partition may change during the evolu- 
tion of the UFO from its origin near the accretion disc 
to the scale corresponding to the smallest cell resolution 
in our simulations (2pc). 

A crucial ingredient in these simulations is the two- 
phase ISM, which consists of a warm (T ~ 10 4 K) phase 
and a hot (T ~ 10 T K) phase. The hot phase is in hy- 
drostatic equilibrium in a two- com ponent gravitation po- 
tential, d escribed by the sum of a Hernquist (1990) and 
a NFW QNavarro et al. 1996) density profile with core 
densities and scale heights (n c ,r c ) = (200 cm -3 , 2 kpc) 
and (3 cm -3 , 20 kpc), representing the baryonic and dark 
matter components, respectively. 

The w arm phase ISM i s initi alized using the algo- 
rithm by |Lewis fc Austin| ( |2002[ ) , which generates a 3- 
dimensional random density distribution that simultane- 
ously satisfies single-point lognormal statistics and two- 
point fractal statistics. These statistical properties are 
consistent with those observed of molecular clouds in 
our G alaxy ( jKainulainen et aL 2009; R oman- Duval et ah] 
2010) and simula tions ol sell-r egulated multi-phase ISM 
turbulence ( Wada fc Norman|2001[ ). We choose the same 
statistical parameters lor the fractal cube as those used 
by SB07. 

The clouds in our simulations are initially in pressure 
equilibrium with the hot phase and we investigate two 
cases for the spatial profile of their mean density: one is 
that of an isothermal sphere, the other is that of a tur- 
bulent quasi-Keplerian disc (SB07). The central mean 
densities, n C) o, are 300 cm -3 and 1000 cm -3 , respectively, 
and the total mass of clouds in both cases is ~ 4x 10 8 M . 
Individual clouds have peak densities of a few 10 5 cm -3 
and temperatures < 100 K. The disc is settled in the y 
- z plane and centered at x = (0, 0, 0) with scale height 
h c determined by the combination of the mean sound 
speed, a « lOkms -1 , and the turbulent velocity disper- 



sion, aturb = 40kms -1 , h c = \J (a 2 + cr t 2 urb )/47rG/in C)0 ~ 
40 pc (SB07). Turbulent support is formally introduced 
to ensure a fairly large disc scale height while keeping 
cloud temperatures below 10 4 K, so that the presence of 
sufficiently large, massive clouds within a scale height 
results in strong interactions between the wind and the 
ISM. 

Advected scalar variables distinguish UFO material 
and warm phase gas from each other and from the hot 
phase background. We include non-equ ilibrium, opti- 
cally thin at omic cooling for T > 10 4 K ( [Sutherland fc| 
Dopita|1993 ).To emulate heating due to photoionization, 



A = for T < 10 4 . Cloud cores initially have tempera- 
tures ~ 100 K with densities up to ~ 10 5 cm -3 . Thermal 
conduction, photo-evaporation, self-gravity, and mag- 
netic fields are not included. 

The cubical simulation domain with physical dimen- 
sions of 1 kpc 3 consists of 352 x 512 x 512 cells in a carte- 
sian grid providing a spatial resolution of approximately 
2pc per cell. Due to the nearly adiabatic nature of the 
flow, a restricted one parameter scaling of physical di- 
mensions by up to a factor of 5 is possible (SB07). 

3. RESULTS 



The evolution of the UFO is shown in Figures [TJ[2j and 

tbr three simulations with different initial ISM. Figure [l] 
ows midplane density slices of the expansion of a UF(J 
into a hot atmosphere devoid of clouds, which we shall 
call run A. Figures [2j [3] show simulations of the UFO 
interacting with a two-phase ISM in which the clouds 
are distributed spherically or in a disc, referred to as run 
B and run C, respectively. 

In run A, the UFO expands in a self-similar fashion into 
a single-phase, smooth, hot hydrostatic medium, giving 
rise to a relatively well-defined two-shock structure and 



unstable contact discontinuity surface (c.f. Weaver et al. 
1977| . Confined by the shock- heated ambient medium 



the turbulent flow beyond the reverse shock circulates 
back outside the freely expanding wind toward the plane 
of the galaxy. 

In the common, very early stage of evolution for runs 
B and C, the freely expanding wind interacts strongly 
with the first clouds in its path and is isotropically di- 
verted into sub-streams. Consequently, any dependence 
on opening angle disappears at this point. Within lOkyr 
after the start of the UFO, however, the evolution be- 
gins to differ between the cases for bulge-like and disc- 
like cloud distributions. In the former, the UFO streams 
continue to branch out isotropically and inflate a quasi- 
spherical energy bubble. The flow entirely engulfs and 
ablates the clouds, primarily driving them radially out- 
ward in long, cometary filaments. The diffuse warm fila- 
ments reach speeds of 1000 km s -1 while the colder cores 
are accelerated up to 200 km s -1 . The results from run 



(2005) 



B are similar to those reported by Saxton et al 
and WBU12 for AGN jet feedback. 

Because of the higher densities and filling factors along 
the galactic plane in run C, the secondary UFO streams 
do not disperse the clouds at large disc radii. As the UFO 
breaks out through the center of the disc carrying with it 
an appreciable mass of dense cloud material, it inflates an 
orbed energy bubble above the disc, which sweeps back 
down over the clouds in the outer regions of the disc. 
These clouds are compressed and accelerated toward the 
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Fig. 1. — Midplane density slices of the evolution of a 10 44 ergs _1 UFO in an ISM devoid of clouds (Run A). See the electronic edition 
of the Journal for a color version of this figure. 
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Fig. 2. — Same as Fig^ but for a two-phase ISM with spherically distributed clouds (Run B). 
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Fig. 3. — Same as Fig^ but for a two-phase ISM with clouds distributed in a quasi-Keplerian dis (Run C). 



galactic plane and central BH by the turbulent, ram- 
pressure dominated back- flow in the bubble. The results 
of thi s run are simi l ar to t hose of the simulations by SB07 
and Gaibler et al. (2012) for AGN jets interacting with 
a dense galactic disc. 

Runs B and C demonstrate that the feedback on the 
warm phase of the ISM depends strongly on the spa- 
tial distribution of clouds. In all runs, however, the 
UFO-blown bubble remains in the energy-driven regime, 
despite radiative cooling in the clouds. This is consis- 
tent with the predictions of r ecent analytic models by 
jFaucher-Giguere fc Quataertj ( |2Q12| ), which also justify 



our neglect of inverse- Compton cooling. 

In the following, we use the four quantities to measure 
the efficiency of feedback by the UFO: the mean radial 
velocity, the velocity dispersion, the mechanical advan- 
tage, and the kinetic energy of the clouds. 

We define the density- weighted mean radial outflow ve- 
l ocity of the warm phase, (v TjW ) = Yl </Ww v ■ r/ Y 4> w Pw 
( Wagner fc Bicknell|2011[ ). The evolution of this quantity 
and its outward only (positive) component are plotted 



together with the total velocity dispersion, a to t^ and the 
(45°) line-of-sight velocity dispersion, cri os ,45, as a func- 
tion of time in Fig. p|a). We see that for the case of a 
bulge-like cloud distribution (run B), the velocities of the 
warm phase reach several 100 km s -1 , and keep increas- 
ing for the duration of the simulation. At late stages of 
the evolution, the clouds are predominantly accelerated 
outward ((v r>w ) « (w r ,w,out))> although their radial speed 
never quite reaches the escape velocity of this system, 
which is ~ 450 km s -1 at 0.5 kpc. The velocity disper- 
sions, however, reach values beyond those predicted by 
the M-a relati on, which for the simulated galaxy using 

the relations by Graham (2012) and a black hole mass of 

T 



6 x 10 7 M o is ~ 170 km s l . The values of (v Y , w ) and a 
are comparable to those found in analogous simulations 
of AGN jet feedback (c.f. WBU12). 

In run C, the feedback in terms of radially outward di- 
rected cloud acceleration and cloud velocity dispersions 
is noticeably less efficient. The radial outflow velocity 
peaks early (after 50 kyr) as bulk cloud material is pushed 
out of the galactic disc and then drops throughout the 
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Fig. 4. — Evolution of various quantities that gauge the feedback 
efficiency for simulations with a bulge-like or a disc-like distribu- 
tion of clouds: a) The density- weighted average radial velocity, the 
outward (positive) component of the radial velocity, the total ve- 
locity dispersion, and the line-of-sight velocity dispersion (at 45° 
inclination); b) The mechanical advantage as measured by the to- 
tal or outward-only radial momenta of clouds, c) The warm-phase 
kinetic energy as a fraction of the energy provided by the UFO and 
the ratio of the warm-phase internal energy to kinetic energy. 



rest of the simulation as the energy of the UFO is pri- 
marily channeled into inflating the bubble beyond. In- 
fall brought about by the surrounding overpressure and 
turbulent backflow dominates the dense gas motions af- 
ter 200 kyr resulting in net accretion. The velocity dis- 
persions also saturate well below 150 km s -1 and do not 
reach the value predicted by the M-a relation for this 
galaxy. 

For a given kinetic power, the ratio of the mass outflow 
rate of the UFO to that of the jet is Mtjfo /Mjet ~ 2(r - 
1)//3ufo ~ 1000, where T is the jet Lorentz factor and 
/3 = v/c. The momentum delivered by the UFO is con- 
sequently larger by a factor ~ 2(r — l)/r/3j e t/3uFO ^ 50. 
In an expanding wind, the momentum transfer leading 
to the acceleration of embedded clouds in all directions 
is provided by the sum of the ram pressure and thermal 
pressure integrated over the surface of the clouds. Be- 
cause the surface area increases over time, this system 
exhibits a mechanical advantage greater than unity, and 
care is required when assessing momentum budgets: the 
net (scalar, not vector) momentum may be larger than 



that injected by the UFO over a given time. Panel b) in 
Fig. [4] shows that the mechanical advantage with respect 
to the clouds, defined here as the ratio of the total ra- 
dial momentum of the warm-phase at a given time to the 
total UFO momentum injected up to that time, in both 
runs B and C is much greater than unity, indicating effi- 
cient momentum transfer^] The efficiency is higher than 
that for AGN jets by almost an order of magnitude (cf. 
WBU12). 

Figure pic) shows the evolution in time of the fraction 
of warm-phase kinetic energy to the integrated energy 
injected by the UFO up to time t, i£kin,w/AjFO x ^ an d 
the ratio of the warm-phase internal energy to its kinetic 
energy, ^i n t,w/^kin,w In both runs B and C, the fraction 
of energy imparted by the jet to the warm phase peaks 
early at > 50% and subsequently declines slowly. Over- 
all, however, the energy transfer efficiency, both in terms 
of heating and accelerating the warm phase is higher in 
the case of spherically distributed clouds, compared to 
the case of clouds distributed in a disc. Due to a higher 
mechanical advantage in the first ~ 100 kyr, the energy 
transfer rate is somewhat higher for UFOs than for AGN 
jets. 

4. DISCUSSION AND CONCLUSIONS 

The simulations presented in this work confirm that 
powerful UFOs are capable of generating strong, galaxy- 
wide feedback. Energy and momentum transfer is 
achieved by fast, mass-entrained flows through the 
porous channels of the two-phase ISM, which carry high 
ram pressure to clouds at all locations of the galaxy, 
even in the plane of the disc. In the case of spheri- 
cally distributed clouds, the feedback results in strongly 
heated and dispersed clouds, accelerated outward from 
the galaxy bulge. In the case of clouds distributed in 
a disc, the feedback results in a rapid lift-up of clouds 
from the plane of the disc followed by compression and 
net inflow of warm disc material toward the center of the 
galaxy. Because of the strong interaction of the freely 
expanding wind with the first obstructing clouds, the re- 
sults do not depend on the opening angle of the UFO. 

Only two ISM distributions have been presented here, 
whereas negative and positive feedback efficiencies de- 
pend on a range of ISM parameters. For example, for 
jet-mediated AGN feedback the feedback efficiency de- 
pends strongly on the mean density and mean size of 
clouds, but only weakly on cloud volume filling factor 
(WBU2012). Feedback by UFOs and AGN jets operate 
alike, and there are several reasons to expect this: 1) 
the injected powers are comparable; 2) the ram pressure 
carried by the fast channel flow is comparable because 
its density is primarily determined by that of the swept 
up hot phase; and 3) there is no dependence on open- 
ing angle. Given that the energy and momentum trans- 
fer mechanisms to the ISM are the same for jet- and 
UFO-driven feedback, it is reasonable to expect that the 
efficiency dependencies on ISM parameters for the two 
scenarios are similar. This work shows in addition that 
the spatial distribution of the clouds (e.g., spherical or 
in a disc) affects the feedback efficiency substantially. 

3 Note that this definition of the mechanical advantag e is some- 
what different to that used by other authors, e.g. Faucher-Giguere 
|& Quataert I (120121). 
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One could now further study the dependence on the 
scale-height of the disc or on the original orientation 
of the outflow. One might expect, for example, that a 
higher degree of misalignment of the UFO with respect 
to the galactic disc normal will lead to stronger negative 
feedback. 

The interactions between an AGN wind or jet and 
the ISM lead to heavily mass-loaded outflows, which is 
compatible with the requirement for ram- pressure domi- 
nated AG N outbursts in the simulations by Gaspari et al.| 
( |2Q12b|a| of the feedback cycle that regulates the ther- 
modynamics of cooling flow clusters. On galaxy and clus- 
ter scales, where mass-loaded jets may have decelerated 



plane by the wind may also become star-forma t ion sites, 
a mechanism proposed recently by Silk et al. (2012) to 



explain Milky Way hypervelocity stars. 

The boundary conditions for the inlet are an extrap- 
olation of the observed UFO parameters, whose origins 
are near the central accretion disc. The evolution of the 
UFO from disc wind on scales ~ 10 -4 - 10 -2 pc to out- 
flow on scales ~ pc is in need of future investigation (see, 
e.g., " 



Wada 2012 



to su b-relativistic velocities (Bicknell 1984 Komissarov 



1994 ), the distinction betw een such slow, massive, wide 



jets (Sternberg et al. |2007| ) and AGN winds may be of 
lesser importance as far as AGN feedback is concerned 



flGaspari etlX|2012b| ) . 

While mass outflow and infall rates can be well de- 
termined, our simulations do not contain all the neces- 
sary physics (e.g. self-gravity, molecule formation and 
cooling, photoionization, etc.) to deduce whether star- 
formation in surviving dense clouds is fully suppressed 
or whether in some cases, most likely within the larger 
cloud complexes in disc galaxies, star- formation may in 
fact be induced, due to the overpress urization of the am- 
bient medium (cf. jGaibler et al.|2012| ). The longevity and 
lateral contraction of filaments drawn out of the galactic 



for radiation-driven outflows). A class 
or component of out flows less powerful than UFO s known 
as wa rm absorbers (|McKernan et al.]|2007| Torresi et al. 
|2012| ) seen in soft X-ray absorption, may help constrain 
the modelling of the wind in this regime. It also remains 
to be seen with simulations of the evolution of the out- 
flow on scales of order tens of kpc whether the filaments 
of dense outward- moving gas may be unbound from the 
galaxy potential and contribute to the enrichment of the 
intergalactic medium. 

A subset of the computations were undertaken on the 
NCI National Facility at the Australian National Uni- 
versity. The work was also supported by iVEC through 
the use of advanced computing resources located at 
iVEC@Murdoch. Ralph Sutherland provided us with 
code to generate the initial warm-phase distribution. 
This work was supported in part by the FIRST project 
based on Grants-in-Aid for Specially Promoted Research 
by MEXT (16002003) and JSPS Grant-in- Aid for Scien- 
tific Research (S) (20224002). 
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